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Lattice QCD is the most reliable non-perturbative method in quantum field theory. In the last few years, some problems crucial to 
high energy experiments have been solved. We review some recent work done by the Chinese lattice community. 
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1 Introduction 

Particle physics and gauge field theories are frontiers of 
fundamental sciences of matter. According to Yang and 
Mills, any basic theory of interacting matter should sat- 
isfy gauge invariance. QCD has been accepted to be 
the most successful gauge theory of strongly interact- 
ing particles. In QGD, hadronic matter is composed 
of quarks, and interactions between them are mediated 
by eight massless gluons generated by the SU(3) gauge 
group. Asymptotic freedom of QGD at short distances 
makes the perturbative calculations of high-energy pro- 
cesses possible. At low energy scales, however, QGD 
has many important properties which can not be stud- 
ied perturbatively, like confinement of quarks and glu- 
ons, vacuum structure, chiral-symmetry breaking, glue- 
ball masses, hadronic spectrum and weak interaction pro- 
cesses, and behaviors of hadronic matters at high tem- 
perature or high density. 

Lattice gauge theory (LGT) has developed into a 
promising first principle non-perturbative approach to 
these phenomena. The basic idea El, as proposed by K. 
Wilson in 1974, is to replace the continuous space and 
time by a discrete grid: 
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Gluons live on links U{x,ii) = e'^^J^; u,x ^^yj, , ^ 
quarks live on sites. The continuum Yang-Mills action 
Sg = J d'^x TrF^i,(x)F^,y(x)/2 is replaced by 



^s = iETi-(Up + Ut-2), 



(1) 



around an elementary plaquette. The continuum quark 
action Sq = J d'^x ■(/'(a;) (7^13^ + m)tlj{x) is replaced by 

Sq — rml>{x)Tp{x) 
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^ ^(a;)[/(x,A:)V(a; + M, (2) 



X k = ±l 



where j^k — ^Ik and the last term is added to avoid 
species doubling. Then all the physical quantities are 
calculable through Monte Garlo (MG) simulations with 
importance sampling: 
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E0[C]-(3) 



Here C stands for a gluonic configuration drawn from 
the Boltzmann distribution. Fermion fields must be inte- 
grated out before the simulations, which leads to O and 
Seff. Since the derivatives are approximated by finite 
differences, the gluonic action in Eq.(|l|) has lattice spac- 
ing error of order a^, and quark action Eq.(H) has error 
of order a. To compare with the real world, the contin- 
uum limit a — > should be eventually taken. On the 
other hand, to keep the physical volume unchanged, 
the number of lattice sites should be very large. There- 
fore, the computational task will then be tremendously 
increased. 

Ghinese physicists have been involved in this field 
since early 80's. Although some very interesting re- 
sults were obtained, most oL them were analytical in- 
vestigations (for review, see 13) , due to limited compu- 
tational facilities. Thanks to the rapid development of 
high performance supercomnuters in Ghina in late 90's 
and Symanzik improvements of LGT, this situation has 
changed. Today Ghinese physicists play an active role in 
the lattice community u. Here we would like to give an 
overview of the relevant developments in recent years. 



where Up is the ordered product of link variables U 
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2 Status 

2.1 Zhongshan University 

The interests of our group cover glueballs, QCD at finite 
density, supersymmetry, new Monte Carlo algorithms, 
and construction of parallel computers. 

(a) Glueball spectrum. The spectroscopy of QCD in 
the pure gauge sector, i.e., the glueball masses attract 
considerable attention. In the quenched approximation, 
these glueball are non qq giuonic bound states formed by 
strong self-interactions of the gluons, and their masses 
vary from about 1.4 GeV to 2.5 GeV. A lot of glueball 
candidates observed in BES experiments such as t(1440), 
/o(1520), e/fj{l720) and ^(2230), produced in the J/i; 
radiative decays are within this range. The difficulty in 
experimental identification of a glueball comes from the 
complexity in determining the quantum numbers J^^ 
of these particles. The lattice QCD prediction for the 
glueball spectrum will help the experimental physicists 
in the their search for glueballs. Concerning the lightest 
glueball 0++ , more accurate MC calculations by the IBM 
group on much larger lattices, higher statistics and bet- 
ter algorithm gave Af(0++) « 1.740 ± 0.071 GeV, where 
the infinite volume extrapolation has been made. Ac- 
tive investigation on the improvement of the lattice tech- 
niques are still being carried out so that the errors due 
to various approximations are reduced and the glueball 
masses are more, accurately estimated. We proposed an 
alternative waycl to extract the glueball masses by solv- 
ing the lattice QCD Schrodinger equation. We obtained 
M(0++) = 1.71 ± 0.05 GeV, in nice agreement with the 
IBM data. The advantage of this method, is its potential 
to extract the wavefunction of a glueball. To reduce the 
lattice spacing errors in Eq.(|l|) andJEfl,(^), we also car- 
ried out the Symanzik improvement Q'Lla, which allows us 
to increase the precision and perform the calculation on 
smaller computers. 

(b) QCD at finite density. This investigation is rel- 
evant for cosmology and neutron star phenomenology. 
When temperature or density is sufficiently high, a new 
state of matter called quark-gluon plasma (QGP) is ex- 
pected to form. The goal of Relativistic Heavy Ion Col- 
lider (RHIC) at BNL and Large Hadron Colhder (LHC) 
at CERN is to create the QGP phase, and replay the evo- 
lution of the universe. Although the standard lattice La- 
grangian Monte Carlo method works very well for QCD 
at finite temperature, it unfortunately breaks down at 
finite density (chemical potential), because of the com- 
plex action. For SU(3) gauge theory with quarks, the 
Lagrangian MC methods always lead to an unphysical 
critical chemical potential /ic = in the chiral limit. The 
Hamiltonian formulation does not have such a problem 
and is therefore a promising alternative. Recently, we 



have developed a_Hamiltonian approach to lattice QCD 
at finite density It avoids the usual problem in the 
Lagrangian Monte Carlo method of either an incorrect 
continuum limit or a premature onset of the transition 
to nonzero quark density as /i is raised. We solved it in 
the case of free quarks and in the strong coupling limit. 
At zero temperature, we calculated the vacuum energy, 
chiral condensate, quark number density and its suscepti- 
bility, as well as mass of the pseudoscalar, vector mesons 
and nucleons. We found that the chiral phase transi- 
tion is of first order, and the critical chemical potential 
(dynamical quark mass at /it = 0). This is 
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consistent with ~ -^^^ /3 (where M^^ is the nucleon 
mass at /i = 0). 

(c) Dynamical quarks. Most conventional full QCD 
algorithms are expensive and do not work in the chi- 
ral limit. The microcanonical fermionic average (MFA) 
method works not only on smaller computers, but also in 
the chiral limit, which is very useful for studying sponta- 
neous chiral symmetry breaking and chiral phasa-transi- 
tion. We generahzed the MFA method to QCDEHI. 

(d) Monte Carlo Hamiltonian. In Lagrangian for- 
mulation, only the lowest-lying state can be extracted 
from the correlation function. It is extremely difficult 
to compute the excited states. We proposed a diflpi- 
ent approach: the Monte Carlo Hamiltonian methodlHl, 
designed to overcome the difficulties of the conventional 
approach. TJip method has been well tested in quantum 



mechanics liiJiij using a regular basis in Hilbert space. To 
apply the method to many body systems and quantum 
field theory, stochastic basis has to be used. 

(e) SUSY. Supersymmetry is a promising fundamental 
theory of elementary particles that describes a mapping 
between bosons and fermions. It underlies many modern 
theories of strings and quantum gravity. Lattice inves- 
tigations of supersymmetry may provide some new non- 
perturbative physics beyond the standard model. Some 
numerical simulations have been donell^. 

(f) Parallel computation. The tremendous advance 
in computer technology in the past decade has made it 
possible to achieve the performance of a supercomputer 
on a vpy small budget. We have built a multi-CPU 
cluster EJ of Pentium PCs capable of high performance 
parallel computations at a very good price/performance 
rate (US$7/MFlops). We believe this is the first such a 
machine constructed in academic institutions in China. 
QCD code with an improved action has been imple- 
mented on the machine. 



2.2 Beijing University 

Although Kogut-Susskind fermions and Wilson fermions 
have been extensively used in numerical simulations, they 
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are not free of problems. There has been evidence show- 
ing that those two approaches may give the topolog- 
ical charge or anomaly incorrectly on a finite lattice. 
Kaplan's domain wall fermions and Neuberger's overlap 
fermion formulation have attracted much attention, be- 
cause they satisfy the Ginsparg- Wilson relations, they 
also produce the correct chiral modes, anomaly and topo- 
logical charge. For domain wall fermions there is an extra 
dimension and the lattice size in this dimension has to 
be very large. Thus algorithms suitable for those new 
fermion approaches need to be developed. Liu did some 
interesting workll3 in this area. 

2.3 Institute of High Energy Physics 

In quantum theory, instantons are known to be a classical 
solution to the 9 vacuum. They might play an important 
role in confinement and chiral symmetry breaking. Chen, 
Wu and He calculated the gLueball masses-using classical 
SU(2) gauge configurations E£l and found liH that the in- 
stantons might also give important contributions to the 
glueball spectrum. Wu is well known for his original work 
with Hamber on improved fcrmionic actions in early 80 's. 

2.4 Institute of Theoretical Physics 

Ma investigated the behavior of the gluon propagator 
in lattice QCD with an improved gauge action E3. He 
also used an improved fermionic action to calculate the 
spectrum of light and heavy hadrons. 



2.5 Zhejiang University 

Ji, Li, Ying, and Zhang began numerical simulations of 
LGT in early 80's. Due to limited computing facilities, 
they changed their interest to statistical physics t3. The 
success of Symanzik improvement retrieves their interest 
in LGT. They recently jiid glueball spectrum calculations 
on a smaller computerc3. 



2.6 Beijing, Nankai and Sichuan Universities 

Chen, and Zheng and Zhu are the leading experts on 
cumulant expansiDn_^Q£ LGT and continue their efforts 
in this direction cdOa. Recently, Ren, Zhu and Chen 
successfully used symbolic language to do higher order 
cumulant expansioncJ systematically. 



3 Outlooks 

These years have seen increasing contributions by the 
Chinese lattice community. We believe they will play 
more active role in the near future. 
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